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We show that in a finite one-dimensional (ID) system with diffusive thermal transport described 
by the Fourier's law, negative differential thermal conductance (NDTC) cannot occur when the 
temperature at one end is fixed. We demonstrate that NDTC in this case requires the presence of 
junction(s) with temperature dependent thermal contact resistance (TCR). We derive a necessary 
and sufficient condition for the existence of NDTC in terms of the properties of the TCR for systems 
with a single junction. We show that under certain circumstances we even could have infinite 
(negative or positive) differential thermal conductance in the presence of the TCR. Our predictions 
provide theoretical basis for constructing NDTC-based devices, such as thermal amplifiers, oscillators 
and logic devices. 



In recent years, nonlinear thermal transport, partic- 
ularly in low dimensional systems, is of significant in- 
terest from both fundamental and practical perspectives 
[HE]- For example, thermal rectification has been exper- 
imentally and theoretically studied in many nanostruc- 
tures [3"HTT] and heterogeneous bulk materials [P2"lffl] . 
Negative differential thermal conductance (NDTC), an 
unusual thermal transport phenomenon where the heat 
current across a thermal conductor decreases when the 
temperature bias increases, is an essential element for 
the construction of thermal transistors [T5] and thermal 
logic |16j . and is shown to exist in many non-linear one- 
dimensional (ID) systems [51 [TT1 H31 HTHIM] and vacuum 
gaps |25j . Many mechanisms such as nonlinear interac- 
tions , molecular anharmonicity [XT] [2T] [57] , inter- 
play between the thermal driving force and the thermal 
(boundary) conductance [T7J (T5J HOI 123] > thermal inter- 
faces [23 HZ1 IH] and others [S3] have been proposed to 
explain the existence of NDTC. Interestingly, several nu- 
merical studies [I8H20] [23] have suggested that NDTC 
may vanish as the system length becomes large (ap- 
proaching diffusive thermal transport). However, it has 
not been definitely answered whether NDTC universally 
vanishes for diffusive thermal transport. Besides, the role 
played by thermal interfaces in NDTC has not been well 
studied. Here, we provide a generic and analytic study of 
these issues in ID diffusive thermal transport described 
by the Fourier's law. We prove that NDTC cannot exist 
when the temperature at one end is fixed. However, we 
show that NDTC in this case is still possible if a junction 
with temperature dependent thermal contact resistance 
(TCR) is introduced. Unlike previous theories and sim- 
ulations HQ E3 ITTHM] that dealt with specific toy 
models which are often difficult to access experimentally, 
our predictions provide a generic way towards building 
NDTC-based devices. 
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We consider a general ID system in the diffusive ther- 
mal transport regime whose thermal conductivity k(x,T) 
is a function of the coordinate x and the local tempera- 
ture T(x). The position dependence of the thermal con- 
ductivity k(x, T) is explicitly expressed, since the system 
we consider can have a spatial dependence of structure 
or composition (e.g., strain or mass gradient). This phe- 
nomenological description is valid as long as the mean 
free path (MFP) of heat carriers is much smaller than 
the size of the system, where the microscopic details are 
unimportant. This approach generates analytic results 
regarding the existence of NDTC, and it is instructive in 
system design to pursue the applications of NDTC. 

For a finite ID system which lies in the coordinate 
range [x L ,x R ] (Fig. nj), the local heat current q(x) can 
be calculated from the Fourier's law: 

q(x) = -k(x,T) — . (1) 

For thermal transport without heat sources or sinks, the 
heat current is conserved and the steady state thermal 
transport equation reads 

Once the temperature at two ends of the system are 
given, i.e., 

T(x L ( R A =T l ( r \ (3) 

the temperature profile T(x) is uniquely [25] determined 
by Eq. ^ and the boundary conditions pj, and the re- 
sulting heat current q (independent of x) flowing in the 
system can be computed from Eq. ([!]). 

By applying an infinitesimal variation ST L ( RS> of the 
boundary temperature at one end, i.e., T L ( R ^ is varied to 
T L ( R )+5T L( - R ) while the temperature T R(L 1 at the other 
end is fixed, the resulting temperature profile is varied to 
T(x) + 8T{x). This temperature profile variation 8T(x) 
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(a) No abrupt junction 
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k 2 {x,T 2 {x)) T" 

(b) Single abrupt junction at x J 



(c) Multiple junctions 
Kn\ K2{x,T 2 {x)) 



X L X J x 3... x n 1 X™ X R + X* — X J X 

FIG. 1. Schematics of ID systems without any junctions (a) 
and with a single junction (b) and multiple junctions (c). The 
junctions are indicated by vertical black lines. 
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can induce a variation Sq of the heat current. We define 
the differential thermal conductance (DTC) as 



G 



Sq 



8{T L - T R ) ' 



(4) 



or specifically 



G = G L = ~y , when T R is fixed; or 
ST L 

G = G R = -^g, when T L is fixed. 



(5) 



In the following we consider the cases with and without 
the junctions to discuss the existence of NDTC. 

Systems without abrupt junctions, as shown in Fig. [lja) . 
In this case, if T R ^ is fixed (without loss of general- 
ity, we can assume T L > T R ), we will show that there 
is no NDTC and the heat current q will increase as 
the temperature bias increases by increasing (lowering) 
T L (R) (NDTC could still exist when the temperatures at 
two ends vary simultaneously |29j . however, we limit our 
study in the cases that the temperature at one end is 
fixed) . 

Qualitatively, we first point out that the non-existence 
of NDTC here is a direct consequence of the uniqueness of 
the solution of Eq. ([2| , as we graphically demonstrate in 
Fig.[2|a). Take the case that T R is fixed as an example. If 
the temperature T L is increased to T' L (> T L ), the tem- 
perature profile T{x) (black line) with T{x L ^) = T L ^ 
and heat current q are changed to T'(x) (red line) with 
T'(x L ) = T' L and T'(x R ) = T R and heat current q' . 
First of all, we must have q' ^ q. Otherwise, the first 
order differential equation ^ = — T j about T with 

initial condition T(x R ) = T R would have non-unique so- 
lutions (T(x) and T'(x)), which is not allowed. Secondly, 
q' cannot be smaller than q (proportional to the slope of 
T(x) at x R ), because otherwise there will be an intersec- 
tion of T'(x) (represented by the dashed line) and T(x) 
at some x 1 < x R . We then must have T'(x) = T(x) in 
the coordinate range [x 1 , x R ] due to the uniqueness of 
the solution to Eq. and thus q' = q (contradiction). 
Therefore, we must have q' > q, i.e., the heat current 
monotonically increases with temperature T L when T R 




FIG. 2. A schematic example of temperature profiles of sys- 
tems (a) without any junctions and (b) with a single junction 
at x J when the temperature at x R is fixed and the tempera- 
ture at x L is increased from T L to T' L . The dotted lines in 
(b) would give rise to NDTC. 



is fixed and there is no NDTC. Similar arguments apply 
to the case when T L is fixed. 

We have derived the analytical expressions for the 
DTCs as (Appendix A) 



G L = ,r\ G R = F(x R ),r\ 



(6) 



where 



F(x) = exp 



J 



1 



dn dT 



x l k(x',T(x')) dTdx' 
F(x>) 



dx' 



(7) 



x l k(x',T(x')) 



dx'. 



Such expressions are useful to calculate the magnitude 
of DTCs from the temperature profile without the needs 
to know directly the heat current and its variation (as 
in Eq. Q). They also directly prove the non-existence 
of NDTC here: since F(x) and k are positive, we have 
G L > and G R > 0. 

Systems with a single abrupt junction, as shown in 
Fig. [ljb). We assume that the abrupt junction is lo- 
cated at x J € (x L , x R ). The system can be considered as 
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two subsystems without any abrupt junctions coupled to- 
gether at x J . We suppose that the subsystems in [x L , x J ] 
and [x^x^] have thermal conductivity «i(x, Ti) with 
temperature profile Ti(x) and k 2 (x,T 2 ) with tempera- 
ture profile T 2 (x) respectively. We denote 



T? = T 1 (x J ), Ti = T 2 (x J ), 
T L = Ti(x L ), T R = T 2 (x R ) 



(8) 



At the junction, the two subsystems are coupled through 
a thermal contact resistance (TCR) R J , defined such that 
the heat current q flowing through the system satisfies 



t j 



T 2 J = qR J (T 1 J ,T 2 J ). 



(9) 



We have generally assumed that the TCR depend on two 
temperatures T( and T/. The TCR R J (Tf,T 2 J ) pro- 
vides the complete characterizations of the junction at 
the phenomenological level. 

The possibility of NDTC here can also be interpreted 
graphically. Again, take the case that T R is fixed as 
an example. If the temperature T L (> T R ) is in- 
creased to T' L (> T L ), the temperature profile T(x) with 
T(x L(fi 0) _ rpL(R)^ j unc tion temperatures T( 2 and heat 

current q are changed to T"(x) with T'(x L ) = T' L and 
T'(x R ) = T R , junction temperatures T[ J 2 and heat cur- 
rent q' , as illustrated in Fig. [2jb) . Because of the discon- 
tinuous jump of temperature at x J , we could have q' > q 
if T'i ] > T 2 J (red dashed line) or q' < q if < T 2 
(red dotted line) where the latter situation gives rise 
to NDTC. We derive the conditions of the existence of 
NDTC in more detail below. 

Analytically, the DTCs are (Appendix B) 



G 
with 



L _( dR J \Fi(x L ) G H_( 1+ dR J 



F 2 (x R ) 



&T() R D 



, (10) 



Fi x =exp / j-^^dx' ,i = l,2, 11 

\J x j m (x', Ti (x') ) dTi dx' J 



defined on [x , x ] and [x , x ] respectively, and 



D 



R 



where 



R J + 1 - q 



OR 1 



Ji 



Fi(x') 



L ki(x',Ti(x')) 



dx', J 2 



dR J 
'dfj 



K2 (x',T 2 (x')) 



■h (12) 



dx'. (13) 



If the TCR is independent of the junction temp erat ures 
T] 1 and T 2 , the partial derivatives of R J in Eqs. ( 10 ) and 
([12]) vanishes, and since Fi in Eq. ( 11 ) and J,; in Eq. ( 13 ) 



are positive, thus G L ( R ^ > and there is no NDTC 
Therefore, a temperature dependent TCR is necessary 
for the existence of NDTC. However, as we will see, it is 
not a sufficient condition. 



In the presence of the temperature dependence of 
TCR, we pick a T B such that T B X > \R- x dR J /dTl 2 \, 
where R = (T L — T R )/q is the thermal resistance of the 
whole system including the TCR. Inside the regime de- 
fined by \T L ~ T R \ < T B in the (T L ,T R ) quarter plane 
we have \qdR J /dTf 2 \ < 1 and subsequently R D > and 
G L{R) > . nQ NDTC is displayed in this low bias regime. 
Therefore, a temperature bias exceeding Tb is required 
to observe NDTC, confirming that NDTC is a nonlinear 
thermal transport phenomenon. 

As the temperature bias (\T L — T R \) increases beyond 
T Bl G L(it) could possibly be negative, leading to NDTC. 
We denote the dimensionless quantities 



such that 



G 



X E 
Y E 

X 



dR J \ Ji 8q 



dR J \ J2 
' dT 2 J ) R J 



dq 



(14) 



■h 



Fi(x L ) 



1+X+Y Ji 



Y 



F 2 (x R ) 



1+X + Y J 2 



(15) 



Now there exists NDTC if and only if at least one of 
X and Y is negative, which means that at least one of 
dqjdT( and —dq/dT 2 is negative [3T]. We refer to such 
junctions as those with intrinsic junction NDTC which 
is now necessary and sufficient for NDTC to occur. Thus, 
the existence of NDTC in systems with a single abrupt 
junction is uniquely determined by the properties of the 
TCRs, regardless of the properties of the system away 
from the junction. 

Furthermore, we can formulate the existence of NDTC 
on the X-Y plane: find out the points on the plane that 
correspond to negative G L or G R . NDTC exists inside 
the shaded areas (A, B and C in Fig. [3|, not including 
the boundaries labelled by the thick solid black and red 
dashed lines. Note that we have G L < and G R = 
on the thin dotted line (-1 < X < 0, Y = 0) while 
G L = and G R < on the thin dash-dotted line (X = 0, 
— 1 < Y < 0). We have R D = and thus infinite (±00) 
DTCs on the thick red dashed line [32] . For the points in 
the shaded areas and close to the thick red dashed line, 
we can have very large magnitude of NDTC, useful to 
design sensitive detectors for temperature fluctuations. 

The blue and cyan shaded areas A and B in Fig. [3] are 
not bounded on the X-Y plane. They correspond to one 
of G L and G R is negative and the other is positive, i.e., 
G L G R < which is equivalent to 



XY < 0. 



(16) 



Eq. ( 16 ) includes the situation that G L ^ is infinite (on 
the thick solid red lines inside the second and the fourth 
quarters of the X-Y plane). We can write Eq. (Jl6| in a 
more transparent way by rewriting the temperature de- 
pendence of the TCR R J (T/,T 2 J ) = R J (f J ,f J ) where 
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FIG. 3. The phase diagram on the X — Y plane that present 
NDTC. 



f J = T{ 



Ti and T J 



dR J 




dR J 


q — s 1 

8T J 


< 


q gfj 



<j 9T j 



(17) 



> 0, i.e., the TCR must 



Eqs. (17) implies that 

be dependent on T J . This is physically significant, be- 
cause in thermal transport we have a natural temperature 
ground of absolute zero temperature. This demonstrates 
the drastic difference between the thermal and electrical 
transport, where in the latter case the junction behaviour 
only depends on the voltage difference (not the average 
voltage) across the junction. 

The red shaded area C in Fig. [3] is bounded and it 
corresponds to the case that both 1 — qdR J /dTf and 
l+qdR J jdTi are negative and R D is positive, equivalent 
to 
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OR' 
df J 



1 > 



dR J 



dT J 



9 < 



R J + Ji + J 2 



8R' 



Jl + 



OR 1 



97V 



(18) 



Eq. (|17J and Eq. ((18) imply that 
dR J 



\q\ > ^max 



and 



\q\ > mm 



dR J 



dR J 



dTi 



dR J 



(19) 



(20) 



respectively, suggesting a positive minimum on the heat 
current for the existence of NDTC. This again indicates 



that the existence of NDTC is in the nonlinear regime, 
beyond low heat current. 

At the onset of NDTC, we have either G L = or G R = 
and the magnitude of heat current (vs. temperature 
bias) reaches its local maximum. Taking the case of G L = 
as an example, we have the following variation rates at 
the vicinity of G L = (Appendix B): 



6 J^ f - r il^F 1 ( X L )-G L (J 1 + J 2 )^F 1 (x L ) > 0, 

5(Tj-T«) L 
6fi ~ G h - 



(21) 



When the system enters the NDTC regime of G L < 0, 
the temperature increase ST L of T L is distributed over 
[x L ,x R ] in such a way that the temperature drop over 
[x J ,x R ] is decreasing with increasing T L (which is the 
manifestation of NDTC) while the temperature drop over 
the junction is increasing with T L , as shown in Fig. [2] 

Systems with multiple abrupt junctions can exhibit both 
NDTC and infinite DTCs, but the detailed conditions 
for their occurrence are more complicated. It can be 
proved that NDTC still requires that at least one of the 
junctions possess intrinsic junction NDTC (Appendix C). 
Nevertheless, these junctions can be grouped into a sin- 
gle effective junction with its properties determined by 
the way the junctions are organized (e.g., the order and 
the connection materials) and by the properties of those 
individual junctions, as shown in Fig. [l|c). After iden- 
tifying the effective TCR R^ s , we can treat the system 
with multiple junctions as one with a single junction. The 
discussions in the previous section can be readily applied 
by simply replacing R J with R^ s . 

This procedure also provides us a routine to engineer 
the TCR. For example, we can construct a system com- 
posed of three segments in [x L ,x{], [x^x^] and [x~,x R ]. 
Suppose that [x L ,x{] and [x 2 , x R ] contain the same kind 
of uniform material and the material in [x^x^] is also 
uniform but different. We can have a single effective 
junction with its TCR R J cS = q/(Tf - Ti) where Tf 
is the temperature at xf at the side of [x L ,x{] 



1 2 



IS 



the temperature at x^ at the side of [rr^! 2 ^] an d q is 
the heat current flowing across the effective junction. In 
this way, we have a symmetrical effective junction, i.e., 
R J cS (T 1 J ,Ti) = R J cS (T 2 \T i '). 

In conclusion, we have studied the steady state ID 
thermal transport in the diffusive regime without heat 
sources or sinks. The Fourier's law is applied to calcu- 
late the differential thermal conductance. We find that 
NDTC (in the case that the temperature at one end is 
fixed) cannot exist in systems without any abrupt ther- 
mal junctions. However, we could have NDTC if and 
only if a junction with intrinsic junction NDTC is intro- 
duced. Our predictions provide a theoretical foundation 
to experimentally realize NDTC through careful thermal 
contact engineering, though it remains an open question 



5 



to realize a junction with intrinsic junction NDTC. 
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APPENDIX 

A. Systems without abrupt junctions 

To calculate the differential thermal conductance 
(DTC), we start from the variation of q = 

- K (x,T(x))^: 



5q 



-k(x, T)—ST 
ax 



On dT 
dfdx 



5T. 



We define 



U L ^(x) 



ST(x) 



(22) 



(23) 



which according to Eq. ( 22 1 satisfies 

d uH R) , 1 dndT L(R) r, L ^G L ^ 
dx n(x,T)dTdx k(x,T) 



= 0. 

where i] L = 1 and rj R = — 1. At the boundaries we have 
ST(x L(R) ) = 6T L W w hile 6T(x R ^) = 0, (25) 



since T(x L ^) = T L( - R \ Thus from Eq. the bound- 
ary conditions for Eq. ( 24 ) are 



U L (x^) = 1, U L (x H ) = U"{x h ) = 0, U H {x") = 1. (26) 
Eq. pil is an inhomogeneous linear ordinary differential 



Ok dT 



and a 



L{R) r ,L(R) 



equation, and the coefficients ^jj^^ K (x.T) 
are functions of x only, since T(x) is already formally 
solved from Eq. ([2| and boundary conditions Eq. 
The solution to Eq. pil is 



L rUR) {x)= um) ^ L ) rfW&M f* F(x') b , 



F(x) 



where 
F(x 



exp 



F(x) 



1 



lk(x',T(x')) 



Ok dT 



dx' 



x l k(x',T(x')) dTdx' 
By evaluating Eq. (271 at x = x R , we obtain 



(27) 
(28) 



G L ^ = J"V (fl) \u L ^(x L ) - F(x R )U L ( R \x R ) 



(29) 



F(x>) 



l k(x',T(x')) 



dx'. 



From Eq. ( 26 ) and Eq. ( 29 ) , we have 



G L = ,r\ G R = F{x R )J~ 1 . 



(30) 



(31) 



B. Systems with a single junction 

For a system composed of two segments lying in 



c L ,x J ] and [x J ,x H ], we denote the temperature pro- 
files Ti{x) and T 2 {x) and the thermal conductivity 
K\{x,T\[x)) and K2(x, Ti(x)), receptively. If a heat cur- 
rent is flowing in the system, the temperature profiles 
satisfy 



q = — Ki(x, T\[x))—. - — , x < x < x , 
da; 

dT 

q = -K 2 (x,T 2 {x)) — — , x J < x < x R 
dx 

with the boundary conditions 



(32) 



T( - T( = qR J (T(,T 2 J ),T 1 {x L ) = T L ,T 2 {x R ) = T R . 

(33) 

Applying the variation ST L ( R ^ of the boundary temper- 
ature at one end, the resulting temperature profiles Ti 
and T 2 are varied to T x (x) + 8T x {x) and T 2 (x) + ST 2 (x) 
respectively. Of course we have 



6Ti(x L ) = ST L , 5T 2 {x R ) = ST 1 



(34) 



The heat current q is varied to q + Sq. The junction 
temperatures Tf and T 2 are varied to T{ J +ST^ and T/ + 
8T2 respectively. We then define the following functions 



tjL(R) = ST 1 , tjL{R) = ST 2 



(35) 



on [x L ,x J ] and [x J ,x R ], respectively. They satisfy the 
following equations 

r^ R ^ R ) = - Ki dU ^ m -^iu^\i = l,2 (36) 
dx a Li dx 

and boundary conditions 

Ut(x L ) = 1, U R (x L ) = U 2 L (x R ) = 0, U R (x R ) = 1. (37) 
Their solutions are 



Fi{x) 



Fi(x) J x .j Ki(x',Ti(x')) 



where 

F 4 (a:) = cxp 



dx' 
(38) 



I 



dKi dT 



j K i (xr,T i (x'))dT i te' 



dx') , i = 1,2. 



(39) 
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By evaluating Eq. (38) at x L for i = 1 and at x R for 
i = 2, we have 



,/ x j Ki(x',Ti(a:')) 



(40) 



' ,j /c 2 (x',T 2 (a;')) 



From Eq. (|40j), we can also readily write down 

5(t? - r/) 



8T L 

s(t 2 j - r fl ) 



= U 1 L (x J )-U 2 L (x J ) 
= F 1 (x L )-G L (J 1 + J 2 ), (45) 
= C/ 2 L (.x J ) = G i J 2 . 



C. Existence of NDTC for systems with multiple 
junctions 



The variation of the first equation in Eq. ( 33 ) gives 



U^ iR \x J ) - U^ R \x J ) = V L ^G L ^R J 

. dR J tt l{r), jx . dR J TT UR)/ jx 



(41) 



By inserting Eq. (40) into Eq. (41), we finally get the 
DTC 



G L 



dR J \ F^x 1 



dT( ) R D ' 

G«-(l + a dRj \ F2{xR) 
G -\ 1 + 9 ffTj) R° ' 



(42) 



where 



with 



" Fl{X>) dx'J^r F2(X,) dx' (44) 
cL nix',Tix')) ' 2 -J xJ Kix'Mx')) ' 



To prove that the existence of NDTC for systems with 
multiple junctions requires that at least one of the junc- 
tions has intrinsic junction NDTC, we prove the converse 
by induction, i.e., there is no NDTC if none of the junc- 
tions has intrinsic junction NDTC. 

Assume that the system lying in [x J ,x R ] contains an 
arbitrary number of junctions and the DTCs for this 
system are non-negative. We now add a new segment 
(with no junctions within the segment) lying in [a; i ,a;' 7 ] 
to the existing system. We denote T ,i (- R ) as the tem- 
perature at x L ( R \ We assume that the new junction 
at x J with TCR i? J (T 1 J ,T 2 7 ) has no intrinsic junction 
NDTC. Here T( and T 2 J are the temperatures at x J at 
the side of [x L 1 x J ] and [x J ,x R ] respectively. Suppose 
we raise the temperature T L infinitesimally to T L + 8T L 
(ST L > 0) while keeping T R fixed. The temperature 



Tf, 2 ) i s then varied to T^ 2 ^ 



■ 8T^, 2 s and the heat current 



is changed from q to q + Sq. At the junction, we have 
Sq = G{ 6Tf - G^STi where G( {2) > are the junction 

DTCs. On the other hand, the system in [x J ,x' R ] has 
non-negative DTCs, i.e., Sq = G B ST 2 T where G B > 0. 
We thus have G{ 8T^ = {G J 2 + G B )5T 2 J . Of course 
STf cannot be negative. If either G 2 or G B is positive, 
we will have 5T 2 J = G{6T^/(Gi + G B ) > and thus 
Sq = G B 5T 2 > and there is no NDTC. If both G 2 and 
G B are zero, we will have 5q — G'lSTf > and there is 
no NDTC. 
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